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Abstract 



■ The intent of this paper is to describe large scale asymptotic geometry of STIT 
Ph . tessellations in M. d , which form a rather new, rich and flexible class of random tes- 
sellations considered in stochastic geometry. For this purpose, martingale tools are 
combined with second-order formulas proved earlier to establish limit theorems for 
STIT tessellations. More precisely, a Gaussian functional central limit theorem for 
the surface increment processes induced by STIT tessellations relative to an initial 
time moment is shown. As second main result, a central limit theorem for the total 
edge length/facet surface is obtained, with a normal limit distribution in the planar 
case and, most interestingly, with a non-normal limit showing up in all higher space 

■ dimensions - including the practically relevant spatial case. 
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1 Introduction and results 

Random tessellations (or mosaics) of lR d are locally finite families of compact convex ran- 
dom polytopes, which have no common interior points, cover the whole space and form 
a central object studied in stochastic geometry, spatial statistics and related fields, see 
[1, 2, 4, 7, 8, 9, 10] to name just a few. However, there are only very few mathematically 
tractable models and the most prominent examples include hyperplane and Voronoi tessel- 
lations, where most often the Poisson case is considered. A new class, the so-called STIT 
tessellations, was introduced recently in [13, 14, 15, 16] and has quickly attracted consider- 
able interest. These tessellations clearly show the potential to become a new fundamental 
reference model for both, theoretical and practical purposes. Whereas most research on 
random tessellations in the last decades was about mean values and mean value relations 
(see [17] for the recent state of the art), modern stochastic geometry focusses on distribu- 
tional aspects [1, 2, 8, 9, 10] and limit theorems [7]. In contrast to the tessellations studied 
so far, the STIT model has the additional feature of arising as a result of a spatio-temporal 
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dynamic construction. From this point of view, limit theorems for STIT tessellations are 
particularly interesting. As recently pointed out by the first author in [18] , we expect that 
the large scale asymptotic of dynamic models will become of great importance in stochastic 
geometry. 

Let us recall the basic construction of tessellations that arise as a result of repeated cell di- 
vision. To this end, let A be a translation-invariant measure on the space TL of hyperplanes 
in R d . Thus, A factorizes under the usual polar identification, i.e. 

A:=£ + ®K, (1) 

where £ + is the Lebesgue measure on the positive real half-axis and where 1Z is a probability 
measure on the unit sphere S^-i- Throughout this paper we always require that the 
support of 1Z spans the whole space, i.e. that span(supp(7£)) = M. d , and we say in this 
case that A is non-degenerate. Further, let t > be fixed and let W C M. d be a compact 
convex window with interior points in which our construction of a random tessellation 
Y(tA, W) is carried out. In a first step, we assign to the window W a random lifetime. 
Upon expiry of its lifetime, the primordial cell W dies and splits into two polyhedral sub- 
cells W + and W~ separated by a hyperplane hitting W, which is chosen according to 
the normalized distribution A. The resulting new cells W + and W~ are again assigned 
independent random lifetimes and the entire construction continues recursively until the 
deterministic time threshold t is reached (see Figure 1 for an illustration). The described 
process of recursive cell divisions is called the MNW- construction in the sequel and the 
resulting random tessellation constructed inside W is denoted by Y(tA, W), as mentioned 
above. 

In order to ensure the Markov property of the above construction in the continuous-time 
parameter t, we assume from now on that the lifetimes arising in the MNW-construction 
(including the initial window W) are exponentially distributed. Moreover, we assume that 
the parameter of the exponentially distributed lifetimes of individual cells [c] equals A([c]), 
where [c] := {H G TL, H H c ^ 0} stands for the collection of hyperplanes hitting cell 
c. In this special situation, the random tessellations Y(tA,W) fulfil a stochastic stability 
property under the operation of iteration of tessellations and are for this reason called 
random STIT tessellations, see Section 2 below for details. 

Having studied the first- and second-order properties of STIT tessellations in [20, 21], 
we consider in this paper the central limit problem. This problem will be considered in 
two closely related settings, interestingly leading to results of very different qualitative 
nature. First, we shall focus our interest on the residual length/surface increment arising, 
respectively, as cumulative length or surface area of the cell-separating (d — l)-polyhedral 
facets born after a certain fixed time in the course of the MNW-construction recalled 
in Section 2. In this set-up we shall establish a central limit theorem with a Gaussian 
limiting variable. Next, we shall pass to the total length/surface, taking into account 
also the polyhedral segments/facets born at the very initial big bang stages of the MNW- 
construction, as descriptively termed in [13]. It turns out that, whereas in dimension 2 the 
Gaussian convergence, is preserved, this is no more the case for dimensions 3 and higher, 
where non-Gaussian limits arise. This apparently surprising phenomenon is in fact due to 
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the influence of the big bang phase in the MNW-construction itself, which is negligible in 
two dimensions but turns out to be crucial in higher dimensions. As already emphasized 
above, it is interesting to see that in terms of the time incremental MNW-construction of 
STIT tessellations our central limit theorems are more reliant on independences arising in 
construction time than those of spatial nature. 

We are now going to describe some of our limit theorems in detail. For W as above, we 
put Wr := RW for R > and let A be some fixed translation-invariant measure on %. 
Therefore and in order to simplify the notation we will write from now on Y(t, W) instead 
of Y(tA, W). Our first limit theorem deals with the total surface area Vo\d-i(Y(tA, W)) of 
the polyhedral facets constructed by the MNW-construction within the time period [s , 1], 
where s > is some positive initial time moment. 

Theorem 1 For each s > 0, the random variable 

^[(Vowni, w R )) - EVoi d _ 1 (y(i, w R ))) 

-(Voi d _i(y( So , w R )) - EVoi d _ 1 (y( So , w R )))\ 

converges, as R — > oo ; in law to jV(0, VV(VoLj-i, A) s 1 ^ d ds), a normal distribution with 

mean and variance IV(VoLj-i,A) j\s 1 ~ d ds, where VV(VoLj-i, A) is explicitly given by 
(9) or alternatively (17) below. 

This statement cannot be extended to s | 0, as would be of interest as potentially yielding 
a Gaussian limit for the (centred and suitably normalized) total edge length/surface area 
VoLj-i(V(1, Wr)). The problem is that the variance integral VV(VoLj-i, A) J so s 1 ~ d ds di- 
verges at 0. However, this difficulty can be overcome for d = 2 but not for d > 2. Indeed, in 
the planar case the asymptotic behaviour of the total edge length turns out to be Gaussian. 

Theorem 2 We have for the STIT tessellation Y(l) in the plane with 1Z given by the 
uniform distribution on Si, 

[Vol 1 (y(l, W R )) - EVol 1 (y(l, Wr))} =}► A/-(0, vr Vol 2 (W)), 



RV^gR 

where =^> means convergence in law. 

In fact, Theorem 1 and Theorem 2 are direct consequences of our much stronger functional 
central limit theorems, Theorems 4 and 5 below. Moreover, we will see there that the 
Gaussian convergence in Theorem 2 is true even for any non-degenerate directional distri- 
bution 1Z. 

For space dimensions d > 2 we claim that the Gaussian convergence cannot be preserved. 
Even though we are able to show this fact for all W and translation invariant A by establish- 
ing non-Gaussian tail decay, for simplicity and in order to keep the argument transparent 
we only give a proof for an more easily tractable particular case, postponing the study 
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of more involved properties of the resulting random field to a future paper. We take 
W = [0, l] d and consider 



where e^, i — 1, . . . ,d are vectors of the standard orthonormal basis for M. d and S rCi+e ± 
is the unit mass concentrated on the hyperplane orthogonal to in distance r from the 
origin. 



converges, as R — >■ oo ; to a non-Gaussian square-integrable random variable E(W) with 
explicitly known variance given by (23) below. 

The plan of the paper is as follows: In the next section we recall some properties of STIT 
tessellations needed for the proofs of our limit theorems. We also recall there some of 
the facts from [20, 21] in order to keep the paper self-contained and present the exact 
statements of our functional central limit theorems. The proofs of our results are the 
content of Section 3. 

We would like to remark that an extended version is available online [19] and, moreover, 
that the results in the present paper form the basis of our subsequent work [22, 23]. 

2 Background material and statement of the func- 
tional limit theorems 

We start by rephrasing some of the properties of the STIT tessellations Y(tA, W) as defined 
in the introduction, the proofs of which may be found in [16]. 

• Y(tA, W) is consistent in that Y(tA, W) n V = Y(tA, V) for convex V C W and thus 
Y(tA, W) can be extended to a random tessellation Y(tA) in the whole space M. d . 

• Y(tA) is a stationary random tessellation, i.e. stochastically translation invariant. 
If, moreover, A is the unit-density isometry-invariant hyperplane measure A; so , or 
equivalently if 1Z in (1) is the uniform distribution Vd-i on Sd-i, then Y(tAi SO ) is 
even isotropic, i.e. stochastically invariant under rotations around the origin. 

• Y(tA) is stable under the operation of iteration, denoted by EEL This is to say 



For this reason, Y(tA) is called a random STIT tessellation in this case. This prop- 
erty was discussed in detail in [20] and we refer to this paper and the references 
cited therein for further discussion, because our arguments do not explicitly use the 
stochastic stability but its consequences. 




(2) 



Theorem 3 For d > 2, A as in (2) and W = [0, l} d , 

R 2 ^ 1} [Vo\ d ^(l, W R ) - EVold_i(l, W R )\ 



Y(tA) = m(Y((t/m)A) EE! • • • EE! Y((t/m)A)), 



m = 2, 3, . . . . 
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Figure 1: Realizations of a planar and a spatial stationary and isotropic STIT tessellation 
(kindly provided by Joachim Ohser and Claudia Redenbach) 

• The surface density of Y(tA), i.e. the mean surface area of cell boundaries of Y(tA) 
per unit volume equals t. In particular, the mean surface area of facets arising in the 
MNW-construction during time [0,t] within a compact convex W C W 1 with interior 
points is given by tVo\d(W). 

• STIT tessellations have the following scaling property: 

tY(tA) = Y(A), 

i.e. the tessellation Y(tA) of surface intensity t upon rescaling by factor t has the 
same distribution as Y (A), the STIT tessellation with surface intensity 1. 

• STIT tessellations have Poisson typical cells, which is to say, the interior of the typical 
cell TypicalCell(y(tA)) of Y(tA) coincides in distribution with that of a Poisson 
hyperplane tessellation PHT(tA) with intensity measure tA, see the discussion after 
Theorem 1 in [20] or [15]. 

For the nonspecialized reader let us remark that the typical cell of a tessellation is what 
we get when we choose equiprobably a cell of the tessellation at random out of a 'large' 
observation window. The exact definition makes use of Palm calculus for which we refer 
to [17]. Moreover, a Poisson hyperplane tessellation PHT(tA) is the random subdivision 
of R d induced by a Poisson point process on the space of hyperplanes % having intensity 
measure tA. 

The finite volume continuous-time incremental MNW-construction of random STIT tes- 
sellations, as discussed in the introduction, clearly enjoys the Markov property in the 
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continuous-time parameter, whence natural martingales arise, which will be of crucial im- 
portance for our further considerations. In fact, this observation was the starting point 
of [20], where a class of martingales associated to STIT tessellations was constructed. In 
order to streamline our discussion we do not repeat the full theory here, but rephrase the 
martingale property of two stochastic processes, on which the proofs of our limit theorems 
are based on. To this end, let Y be some instant of Y(tA, W) and let <${■) be a measurable 
facet functional of the form 

(KiO-voW/KW)) (3) 

with n(/) standing for the unit normal to / and ( for a bounded measurable function on 
Sd-i- Moreover, denote the collection of cell-separating (d— l)-dimensional facets, usually 
referred to as (d — 1) -dimensional maximal polytopes, arising in subsequent splits in the 
MNW-construction by MaxPolytopes d _ 1 (F), define E^(Y) by 

MY) ■= E 

/eMaxPolytopes d _ 1 (y) 

and A<p(Y) to be 

MY) := / E <PU)KdH) 
J ^ /eCciis(ynH) 

with Cells(Yflif) standing for the set of (d — 1) -dimensional cells of the tessellation Y(lH 
induced by the intersection of Y with hyperplane H. Let us further introduce the bar 
notation S^(Y) for the centred version S^(Y) — EE^(Y) of S^(Y). Then, we have (see 
[20, 21]) 

Proposition 1 The two stochastic process 

^(Y(tA,W)) and Ej(y(tA, WQ) - / A^(Y(sA, W))ds (4) 

Jo 

are both Q t -martingales, where ^s t stands for the filtration generated by (Y(sA, W^)) < s <t- 

In particular, see [11, Thm. 4.2], the martingale E^Y^A, W)) has its predictable quadratic 
variation process (E0(Y(-, W))) absolutely continuous and given by 

(^(Y(;W))) t = f A^(Y(sA,W))ds. (5) 
Jo 

Beside these martingale tools, we will also make use of the following formula for the variance 
Var(E (Y(tA, W))) of E (Y(tA, W)), W C R d compact, convex and with interior points, 
established in full generality in [21] in order to calculate the variance of the limit random 
variable of our non-Gaussian limit theorem. 
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Proposition 2 For any non- degenerate translation-invariant A on H and as in (3), we 
have 

Var(E,(y(tA,W)))=/ ( 2 (n(iJ)) / / 1 ~ ^ff^ ^yA(dg). 

^[W] JwnH JwnH IV \l x y\) (g) 

Let us further recall from [21] that the variance of the total edge length of a stationary 
and isotropic STIT tessellation Y(tA iso , Wr) in the plane behaves asymptotically like 

Var(Voli(y(*A iso , W R )) ~ tt Vo\ 2 (W)R 2 logR, R^oo, (7) 

where W is again a compact convex set and RW stands for W dilated by a factor R > 0. 
Indeed, this can be seen from the general statement from Proposition 2 combined with 
integral-geometric tools. Note that the asymptotic variance expression for the total edge 
length is independent of t. However, for all space dimensions > 2, the surface density t 
enters the asymptotic variance expression as shown in detail in [21]. 

We can now turn to the statement of our functional limit theorems, from which Theorems 
1 and 2 are direct consequences of. 

Theorem 4 For each s > the centred surface increment process 

■= 4u2 [^( y (*> w *)) - ^(Y(s , W R ))] 

n / *e[so,i] 

converges in law, as R — >■ oo, on the space V[s , 1] of right continuous functions with left- 
hand limits (cadlag) on [s , 1], endowed with the usual Skorokhod topology [3, Chap. 3, 
Sec. 14], to a time-changed Wiener process 

1 W vv(M)/,* - 1 -« , d»' 

where W(.) is the standard Wiener process and Vw(<f>, A) is given by (9) or alternatively 
(17) below. In particular, 

= ^72 $+(Y(l, W R )) - ^(Y(s , W R ))] 

converges in law to Af(0, Vw(4>, A) s 1 ~ d ds), a normal distribution with mean and vari- 
ance Vw(<j>, A) s 1 ~ d ds. 

We consider now the functional convergence of the total length process in the planar case. 
Write 

t(s,R) := expQlogi? -log log R}(s - 1)) = R 3 ' 1 (log R) 1 ' 3 
and define the total length process 

Theorem 5 The total length process (C^' w ) S £[o,i] converges in law, as R — > 00, on the 
space X>[0, 1] of cadlag functions on [0, 1] endowed with the usual Skorokhod topology, to 
(a/VV(0, A)V\4) s g[o,i] where, again, W(.) stands for the standard Wiener process. 
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3 Proofs 

After having rephrased some background material on STIT tessellations in the previous 
section, we are now prepared to present the proofs of our limit theorems. Let us briefly 
recall that we will deal with a fixed translation-invariant hyperplane measure A and that 
for this reason we will write Y(t) instead of Y(tA) without confusion. Moreover, we fix 
some compact and convex set W C M. d having interior points and write Wr = RW for W 
dilated with a factor R > 0. Moreover, recall that the face functionals we are dealing with 
have the representation (3), that S^,(F(t)) was defined in (4) and that the bar notation 
S^(y(i)) stands for the centred version S^(y(i)) — EE^(Y(t)). We start now with the 

Proof of Theorem 4. Notice first that, because of AflWa]) = RA([W]), 

l- d A^(Y(l,W R )) = ± I j^e^{H)) VolLMMdH) 

feCeUs(Y(l,W R )nH) 



= I ^C 2 (fi(tf)) £ Vol^/Mdff). (8) 

J \ W ^ feCe\\s(Y(l,W R )nRH) 

We claim that, upon letting R — > oo, this converges in probability to 

f 9,^, u , xEVol^ JTypicalCellfrri) n if)) , . 

v i tun f > E Vo ^- 1 (TypicalCell (K( 1 )n.\ „ 

where -u -1 is the orthogonal complement of u G Sd-\ and 7?. is the directional distribution 
of the stationary STIT tessellations Y(t) as given in (1). To see it, recall that Y(l)nRH is 
a STIT tessellation in RH for each R > and H G "H. Thus, applying [17, (4.6) and Thm. 
4.1.3] to this tessellation and the fact that E Vol rf (TypicalCell(F(l) n -U" 1 )) is the same as 
the inverse cell intensity of the tessellation Y(l) n -U" 1 , see (10.4) ibidem, we have 

lim — 1-i-E V Vol^ (/) 

i^oo Rd-l d—l\J 
feCe\\s(Y(l,W R )nRH) 

VI fw n m E Vog-i(1VpicalCeU(y(l) n g)) 

= vou-x^ n g) 1 Vold _ i(TypicalCell(y(1) n H)) ■ (io) 

Next, we observe that Y(1,W R ) n i?if = Y(l) n i? - H (H n W), where - H is the scalar 
multiplication relative in ii, that is to say H 3 R rr = p/f(0) + i?(x — p/r(0)), x G 
H with standing for the orthogonal projection on H. Thus, applying the recently 
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developed strong mixing and tail triviality theory for STIT tessellations, especially [12, 
Thm. 2], noting that tail trivial stationary processes are ergodic [6, Prop. 14.9] and then 
using the standard multidimensional ergodic theorem, see e.g. Corollary 14. A5 ibidem, to 

W^J2 f eCciHY(i)nR- H (Hnw)) Vol Li(/)> we S et from ( 10 ) that 

lim V Vol' ,(/) 

feCc\\s(Y(l,W R )nRH) 

Vnl (w n m E VolL^TypicalCelKril) n H)) 
d " U J EVol d _ 1 (TypicalCell(F(l) n H)) 

in probability. Putting this together with (8) and integrating over [W] yields 

lim 1^A^(Y(1, W R )) = V w (<f>, A) in probability, (11) 
as required. Note now that by the scaling properties of Y(s, W R ) and 4> 2 for s > we have 



±- d A, 2 (Y(s,W R )) £ w L- I A^(Y(l,W aR )) £ ^ ^^(¥(1, W sR )). 
Thus, combining (11) with the scaling relation (12) we get 



(12) 



lim ±- d A^(Y(s,W R )) = -Lv^^A) (13) 

in probability uniformly in s G [so, 1]. This crucial statement puts us now in context 
of the general martingale limit theory. Indeed, using Propositon 1 we see that <S^ '^ = 
-^72 [S^,(F(1, W R )) — S^(F(s , W R ))] is a martingale with absolutely continuous predictable 
quadratic variation process 

(S% W )t = jf ^A^{Y{s, W R ))ds, (14) 
see [11, Thm. 4.2]. In these terms, (13) yields for each t 

lim (<S*f > t = / -Lv w (<f>,A) in probability. (15) 

R-+00 J SQ S d 

We now want to apply the martingale functional central limit theorem. Whereas this is 
well known for continuous martingales, we need a version for martingales in the Skorokhod 
space T>[sq, 1]. In this paper, we will make use of the version formulated as Theorem 2.1 
in the survey article [24]. In order to apply this theorem, several conditions have to be 
checked. Condition (ii.6) in [24, Thm. 2.1] is just (15), whereas condition (ii.4) there 
is trivially verified, because the predictable quadratic variation (S^ ' w ) has no jumps by 
(14). It remains to check the condition (ii.5) ibidem, which is that the second moment of 
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the maximum jump J(S^ W ; 1) of the process (<5^'^) se [ So ,i] goes to as R — > oo. More 
precisely, 

<j( qR,W . i \ _ \ q R,W q R,W i _ ,. 

^l^so,- SU P I^V* -^s ,t-\, o SOjt _ - limJ> S0 S , 

s <t<l S T* 

and we have to check that EJ 2 (S^ '; 1) — >■ 0, as i? — >■ oo. To this end, note first that, with 
probability one, J{Sf^ w ] 1) is bounded from above by a constant multiple of Rr d l 2 times 
the id — l)-th power of the diameter of the largest cell of F(s , Wr). Since the typical cell of 
Y(so) is the same as that of a Poisson hyperplane tessellation with intensity measure sqA 
(see Theorem 1 in [20] or Section 2 above), we conclude by standard properties of Poisson 
hyperplane tessellations that the expected number of cells in Y(s , W R ) with diameter 
exceeding D is of the order 0(R d exp(-D)), since A has been assumed to have its support 
spanning the whole of R d . Summarizing, setting u = D d ~ l Rr d l 2 , we are led to 

V(J(S%™ ; 1) > u) = 0(R d eM-R d,(2d ~ 2) u 1,(d - 1) )). (16) 
Clearly, (16) is much more than enough to guarantee that 

which gives the required condition (ii.5) of Theorem 2.1 in [24]. Upon a trivial time change, 
this theorem yields now the functional convergence in law as stated in our Theorem 4. □ 

Before turning to the proof of Theorem 5 we provide an alternative formula for the factor 
Vw(cj),A). Readers not specialized in convex or stochastic geometry could also skip this 
alternative representation and directly jump to the next paragraph, because Proposition 3 
will not be used in the sequel. However, having such a more explicit variance expression 
could be useful for other purposes and was already used in our subsequent work [23]. We 
denote, as in [20], by II the associated zonoid of a Poisson hyperplane tessellation with 
intensity measure A, by 11° its dual body and by 1Z the directional distribution of the STIT 
tessellation from (1), see [17] for the precise definition of IT. 

Proposition 3 We have 

Vw{4>, A) = Yo\ d (W)^-^ [ ( 2 (u) VoWOlltOWdu), (17) 

where Ulu 1 - stands for the orthogonal projection of U onto the hyperplane u 1 - , and where the 
polar body (IIIm- 1 ) is considered relative to u^. In the isotropic case, i.e. when 1Z = va-i, 
the uniform distribution on the unit sphere Sd-\, this reduces to 

W(0,A iso ) = VoLKW/^-V-lr (^) d V (£f " J s C 2 {u)v d ^{du). 

In particular for ( = 1, W = Bf the unit ball and d = 2 and d = 3 we conclude the exact 
values 

32 

V B 2(Yo\ 1: A iso ) = rr 2 and V Bf (Vo\ 2 , A iso ) = — tt 2 . 
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Proof of Proposition 3. At first, [5, Cor. 3.7] provides a general formula for the 
second moment of the volume of the typical Poisson cell of a stationary Poisson hyperplane 
tessellation in M d . In terms of the zonoid II it reads 



,2,-,,. ■ nmUT , . n » d\ Voi d (n c 



E Vol 2 d (TypicalCell (PHT ( A) ) ) 



2 d Vold(n) ' 



where we have used formula [17, (4.63)]. Moreover, the mean volume of TypicalCell(PHT(A)) 
is given by 

EVol d (TypicalCell(PHT(A))) = — 1— 

Void (II) 

according to [17, Thm. 10.3.3 and (10.4)]. Using now Eq. (4.61) ibidem and the fact 
that STIT tessellations have Poisson typical cells and replacing d by d — 1 in the last two 
formulas we obtain (17). The precise value in the stationary and isotropic case can be 
calculated from the fact that in this case, II is a ball with a known radius, see [17]. □ 

We are now going to present the 

Proof of Theorem 5. Note first that 

t(0,R) = ^-, t(1,R) = 1, — r(s,R) = r(s,R)[\ogR-\og\ogR]. (18) 

Thus, defining the auxiliary process 

M*> W = M S := 1 Tj\2*(Y(T( 8 , R), W R )) - E»(y(r(0, R), W R ))} 

Ry/ log R — log log R 

and using (4) with W R := RW and under variable substitution s := t(u,R) and t :— s 
with LHS variables corresponding to the notation of (4) and RHS to that used here, we 
see that, by (18), 



(M s )l =0 and (m 2 s - J* 1^±A^(Y(t(u, R), W R ))du 



«e[o,i] 



are QV( S) #)-martingales. In particular, see [11, Thm. 4.2], the predictable quadratic varia- 
tion process (M) s is given by 

(M) s = £ Ti^l A4>2 (Y(r(u,R),W R ))du, s e [0,1]. (19) 

Repeating the argument leading to (13) we see that 

lim T^JV A(t>2 (Y(T(s,R),W R )) = VW(0,A) (20) 
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in probability, uniformly in s G [0, 1]. Note that the uniformity in s comes, as in the case 
of (13), from the relation (12) implying that, in distribution, all instances of the LHS 
for different values of s are just scaling instances of the same object R~ 2 A ( f,2(Y(l,W^)) 
for R = R/t(s,R) and thus, in terms of the considered convergence in probability to a 
deterministic limit, we are just dealing with a single asymptotic statement. Consequently, 
by (20) and in full analogy to (15), 

lim (M) s = [ V w ((f), A)du — sV w ((j), A) in probability. (21) 

R-+00 J 

Thus, we are again in a position to apply [24, Thm 2.1] yielding the functional convergence 
in law, as R — > oo, in ©[0,1] of (M s ) se [ 01 ] to (^/Vw(4>, A)W s ) se [o,i]. Indeed, condition 
(ii.6) there is just (21), condition (ii.4) is trivial in view of (19), whereas the condition 
(ii.5) is verified by noting that, with probability one, J(M; 1) = r O(R diam( W) ) = 
0(l/Vlog-R), so that in particular lim^oo Ej7" 2 (M; 1) = 0, as required. Denoting now by 
C R ' W the correction term (i? v 45g _ R)- 1 S^(F(r(0, R), W R )) such that 



r n,w r R,w , / log # -log log i? ^ 
C ° = ° + V bgl2 Ms ' 

noting that log R — log log R ~ log R and that, by the scaling property of STIT tessellations 
and by (7), 

Var(C^ R ) = ©([^(log^-^i^/aog^^Oog^^loglogi?)]) = OQoglog R/ log R), 

we see that the processes M s and C R ' W are asymptotically equivalent in V[0, 1] as R — > oo. 
This completes the proof of Theorem 5. □ 

In the context of proof of Theorem 5 it should be remarked that the 'negligible correc- 
tion term' C R ' W has its variance of order 0(loglogi?/ log/2) and thus indeed tending to 0, 
but extremely slowly. Consequently, although the Gaussian CLT holds for C R,W , it is quite 
natural to expect that the convergence rates are extremely slow, conjecturedly logarithmic. 
This is due to the fact that dimension 2 is the largest dimension (critical dimension) where 
the Gaussian limits are still present. In dimensions 3 and higher there is no Gaussian CLT 
and the 'correction term' analogous to C R,W will turn out order-determining rather than 
negligible, as shown by Theorem 3. 



Proof of Theorem 3. Even if in the formulation of the theorem we have used the vol- 
ume functional, we will show the statement in a more general context, where the volume 
VoLj-i(l, Wr) is replaced by a general face functional S^(l, Wr) satisfying (3). 
We claim that the argument from the proof of Theorem 5 cannot be repeated for d > 2. 
Intuitively, this is due to the fact that for d > 2 the variance order of Tl^,{Y{\,Wr)) 
is 0(R 2( - d ~ 1 ^), see below, whereas the variance order of the increment S^,(y(l, Wr))- 
S^(F(so, Wr)), with some time instant < s < t, is 0(R d ) as seen by Theorem 4. 
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Hence, for d > 2 we conclude that even the very first faces born in the MNW-cell-division 
process already bring a non-negligible contribution to the overall variance. Thus, we can- 
not split the whole STIT construction into the warm-up phase (t e [0, R^ 1 log-R] for d = 2) 
with negligible variance contribution and the proper phase unfolding already in a typical 
STIT environment. In fact, the claim is that the CLT does not hold for STIT surface 
functionals in dimension greater than 2! 

Recall that we do not show this fact in full generality for all non-degenerate hyperplane 
measures A and all windows W, but restrict ourself to a particular case, where A is given 
by (2) and where W = [0, l] d . To see the non-Gaussianity, observe first that, by the scaling 
property of STIT tessellations, 

ir^E^ya, W R )) = ^(Y(R, W)), (22) 

which implies that the variance Var(E«^(y(l, Wr))) is of order 0(R 2( - d ~ 1 ^). Indeed, this 
follows directly from the special form (3) of the face functional 4> and the scaling relation 

Y(1,W R ) = RY(R,W). Further, recall that by (4) the process R ^ ^{Y {R,W)) is 
a square-integrable martingale with absolutely continuous predictable quadratic variation 
process given in (5) and, moreover, by Proposition 2 we have 

Va I p*(Y(R,W)))= [ ( 2 (*(H)) / / 1 ~^}~ R ^ [xy])) dxdyA(dH), 

J[W] JwnH JwnH A (NJ) 

which is bounded uniformly in R. Consequently, by the martingale convergence theorem, 
there exists a centered square-integrable random variable such that 

E(W) = lim Y,^Y{R,W)) 

R-^KX) 

a.s. and in L 2 and, moreover, 

Vax(S(W)) = lim Var(E < ^(y( J R, W))) 

= / C^{H)) [ [ —!—dxdyA(dH). (23) 
J[W] JwnH JwnH ^\[ x y\) 

Using now (22) we readily conclude that 

R-^%{Y{l,W)) =>E(W) 
as R — > oo, where =^> stands for convergence in distribution. 

We show now that the variable S(W) cannot be Gaussian. To see it, consider the event 
£n, N > 0, that only hyperplanes orthogonal to e\ have been born during the time [0, 1] 
in the MNW-construction and that their number exceeds N. Observe that, in view of the 
form (2) of A, F(£ N ) = exp(-d) Yl'kLN+i M anc ^ ^ nus 

\og¥(S N ) = -&(N\ogN), (24) 
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where by ©(•) we mean something bounded both from below and above by multiplicities 
of the argument, i.e. ©(•) = O(-) n Q(-) in the usual Landau notation. Further, given the 
fixed collection of all hyperplanes H\, . . . , Hk, k > N, born at times between and 1, on 
the event En, we see that the conditional law of S(W) coincides with that oik — d plus the 
sum of independent copies £1, . . . , Ck+i of E(Wi), . . . , respectively, where Wj, j = 

1, . . . , k + 1, are the parallelepipeds into which W is partitioned by H 1: . . . ,H k . Note that 
the extra k above is the sum of (d — l)-volumes of W n iJj, whereas — d = — EE^(y(l, W)) 
is the centring term. Since Var(£i + . . . + £k+i) — Y2j=l Var E(Wj), which is bounded from 
above by Var(S(VF)) in view of (23), by Chebyshev's inequality we get 

P(6 + • • • + 6 +1 > -2 v / Var(H(^))) > 1 - ^^^f^ > 3/4. 

Thus, in view of (24) 

F(E(W) >N)>lF (^ +2V ^p^y +d ) = exp(-0(iVlogiV)). 

Since Gaussian variables exhibit tail decay of the order exp(— 9(iV 2 )), the random variable 
E(W) cannot be Gaussian, which completes our argument. □ 
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